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SOCIEDADE DE PRÉ-FABRICAÇÃO E OBRAS GERAIS, NOVOBRA, LDA. 
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TEM CANALISAÇÕES DE VAPOR NA SUA FÁBRICA ? 
TEM TRANSPORTE DE GASES? 
TEM TRANSPORTE DE FLUÍDOS CORROSIVOS ? 


Então tem problemas nas tubagens : 


| — Fugas 2 — Corrosão 3 — Calcinação das juntas 
ÃO SEU SERVIÇO — a última conquista da quimica neste sector: 


CROSSLITE — Fita vedante para ligações 


tubulares à base de Polytetrafluoreto de Etilene 


+ Adesão e cobertura perfeita de todas as irregula- 
ridades da rosca, 


+ Insolúvel e absolutamente inerte ao ataque quí- 
mico de qualquer substância. 


+ Não contamina as substâncias em processo através 
da tubagem. 


+ Resiste a temperaturas de — 1000 C a + 3000 €. 

+ Preserva impecavelmente o sistema roscado ao fim 
de meses cu anos de trabalho, 

Um produto de 


HENRY CROSSLEY LTD. — Bolton — ENGLAND 


| bobine: ca. 4 mts, comprimento Todos os acessórios para vedação 
ca. ),5 cm. largura 


Representante SIDROL — Sociedade Importadora de Drogas, Lda. 
EECINSTOO à Calçada de Monchique n.º 3 — PORTO — Telefs. 27826 e 32055 
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- 1ÉvVAISs 405 MELHORES ESTRANGEIROS 
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Quase 40 Anos 
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Bombas centrífugas «PACIFIC» 


Batignolles - Chatillon 


Especialmente destinadas pare: Sob as condições de funcionamento 


— Petróleos seguintes: 

— Indústrias quimicas — Débitos: 1 m*/h a 10.000 mº/h 
— Construções navais — Pressões: altas e baixas 

— Centrais eléctricas — Temperaturas: —-180ºC +600º0 


— Alimentação de caldeiras | 


— Líquidos: corrosivos ou não 
— Etc., etc. q 


Para estudos e orçamentos, consultar o representante em Portugal: 
CIDEX — Comércio Internacional 
e Representações Industriais, Lda. 
Av. Casal Ribeiro, 46-2.º/Telef. 7353 07/LISBOA 1 
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Máxima protecção 
Facilidade de aplicação 
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AGÊNCIA COMERCIAL, L.* 


"* LISBOA EVORA 
T. do Loureiro, 3 R. da República, 93 


Telef. 40054 Telef. 233 63 
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Válvulas de passagem 
Sede de estanqueidade e haste de açu 


: inoxidavel 
Válvulas de cunha, para vapor oxidave 


Placas de obturação articuladas sobre uma 
esfera. Estanqueidade permanente 


2 Valvulas para matérias pastosas 
Sem necessidade de conservação (sem 
prensas-estopa nem tampas de limpeza) 


3 Valvulas de sede obliqua 
para ácidos. Liga cromo-niquel-molibdeno 


4 Valvulas de membrana de borracha 
Sem prensa-estopa 


5 Valvulas de macho 
Com elevada resistência aos ácidos 


6 Valvula de cunha. Sem prensa-estopa 
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Fábricas de Klus (Suiça) 


a frente no progresso da 
construção de válvulas 


Representantes em Portugal. 
Socotel, Lda. 

Rua Sá da Bandeira, 651-4º, Porto 
Tel. 27013 


TECNICA V 


E fio MOTORES INDUSTRIAIS 
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MOTORES PARA TODAS AS APLICAÇÕES INDUSTRIAIS 
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OUTROS MODELOS PARA ENTREGA MUITO RÁPIDA 


MOTORES A 4 TEMPOS - ARREFECIMENTO POR ÁGUA 
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O EMPREGO DO CIMENTO BRANCO... 


permite acabamentos 

mais perfeitos, mais 

duradouros e muito 
mais económicos 


Estude a vantagem do emprego do 


CIMENTO BRANCO LUSO 


Consulte os distribuidores gerais 


No Sul: SCIAL —T. do Corpo Santo, 15 — Telef. 66 61 86 — Lisboa 
No Norte: SCIAL — R. do Bonjardim, 205 — Telef. 25779 — Porto 
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COMPANHIA PORTUGUESA 
DE FORNOS ELÉCTRICOS 


PRODUTORA DE 


FERRO GUSA 
FERRO LIGAS 
FERRO -MANGANEÊS 
SÍLICO-MANGANÉÊS 
FERRO-SILÍCIOS 
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CARBONETO DE CÁLCIO 
CIANAMIDA CÁLCICA 
PASTA PARA ELÉCIRODOS 


E EM MONTAGEM 


GRAFITE 
FÁBRICAS SEDE 
CANAS DE SENHORIM Largo de S. Carlos, 4 2.º — LISBOA 


Telefone 67222 Telefones 325343-355 56-36 89 89 
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Terraços e pavimentos isolados com CELULITE no 
edifício da Comp." de Seguros Fidelidade — P. das Águas Livres — Lisboa 


SETH, LDA. + RUA FILIPE FOLQUE, 10-1.º + TEL. 7301 56 — LISBOA 
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Para impermeabilizar terraços, 
paredes, fundações, etc. 


COLADEM DE 
TACOS DE MADEIRA 


TELEF. 730156  TELEG. EPALDA — LISBOA 


FPABRICA EM SACA TEN 


TECNICA X 


há a 


ATARAAA A 


Mm TA 


TD 


nuunnnnnuaa 
LAS 


MURAL 


ERA 


ETTA 


DA 


E Apesar | Sd 


ND 


Ena | 


= o 


Lataana nu 


| AA, 


£. grande marca italiana 
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Rua Rodrigo da Fonseca, 


Tel. 24818 — PORTO 


LISBOA — Tel. 686072 


TECNICA XI 


Ateliers de Construction Oerlikon—Zurich 50 


Alternadores accionade: por turbinas hidráu- 
licas, para todas as potências e condições de 
exploração, pequenas e grandes velocidades 


Central hidroeléctrica de Birsfelden (Suíça), sobre o Reno, Alternauores de 
26.000 kVA, cada um, 6.600 V, 68 rot/min, 50 Hz. Diâmeiro: 14 metros 
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Ano XXXVII— N.º 331 Outubro de 1963 


C. D. U. 378.962 (064.3) 


CONGRESSO DO ENSINO DE ENGENHARIA 


CONCLUSÕES GERAIS 


Tirar conclusões, bem fundamentadas e que correspondessem à opinião da maioria dos Congres- 
sistas, das comunicações apresentadas e das intervenções registadas durante as sessões de trabalho, 
revelou-se tarefa ingente. É possível que, com mais vagar, se pudesse vir a formular um conjunto 
de conclusões que exprimissem, com mais felicidade, o pensar geral, todavia, sempre à custa 
de um número bastante elevado de conclusões e de dilação grave da sua apresentação a público, 
que já não é pequena. 

Nestas circunstâncias, optou-se por extrair do conjunto das comunicações e, sobretudo, das 
apreciações e discussões que fomentaram, um número limitado de conclusões de carácter geral 
que exprimissem os anseios manifestados e as directrizes apontadas, em vez de conclusões mais 
pormenorizadas mas, certamente, mais controversas. Não é, aliás, de estranhar que tal suceda, 
uma vez que os problemas do ensino têm sido pouco estudados e discutidos entre nós, onde se vive 
muito da tradição, da cópia e do empirismo em tal assunto. 

A Comissão Organizadora do Congresso elaborou um plano de trabalhos, bastante detalhado, 
com 16 temas propostos para estudo, julgando que a indicação de temas objectivos e limitados, 
facilitava a apresentação de comunicações sobre assunto tão complexo como o ensino superior 
de engenharia. O que se verificou na prática mostrou quão inadequado tinha sido aquele modo 
de proceder — raras foram as comunicações que não focaram vários temas de estudo propostos — 
embora não se veja como teria sido possível programar as sessões de trabalho de outro modo. 

Este mal de origem viu-se ainda agravado pelo facto de a discussão ter posto em causa não 
sômente o ensino superior de engenharia, mas todo o sistema de ensino. A complexidade do problema 
ficou bem expressa nas implicações que envolveu e a que não se fugiu, nem seria possível fazê-lo, 
durante os debates que o Congresso provocou. 

Como síntese de tudo o se que passou durante o Congresso, pode afirmar-se ser dever imperioso 
promover urgentemente uma reforma profunda do ensino em todos os escalões, com base em estudo sério, que 
analise não sóas necessidades do País mas também as circunstâncias do ambiente, pois que o problema do ensino 
é, no fundo, um problema humano. 

O Congresso do Ensino de Engenharia formula as seguintes condições : 


I—É de importância vital e muito urgente a remodelação da actual estrutura do ensino 
da engenharia, por não satisfazer às necessidades da vida nacional, implicando tal remo- 
delação a reorganização do ensino em todos os seus aspectos, nomeadamente : 


a) elevação do nível de preparação geral do País por meio da extensão do período de 
escolaridade gratuita e obrigatória, com um imediato alargamento para 6 anos numa 
primeira fase ; 
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b) 


c) 


d) 


f) 


8) 


preparação pré-universitária dos alunos de engenharia com uma índole essencialmente 
formativa e dirigida às cadeiras de base no que respeita à formação técnica dos 
futuros engenheiros (Matemática, Física, Química e Desenho), mas ainda em relação 
à formação integral de futuros responsáveis (Línguas, Filosofia, História, Artes, 
Geografia, etc.) ; 

instituição de graus académicos com um primeiro grau de duração não superior 
a 5 anos; 

reestruturação do ensino em harmonia com a alínea anterior facultando-se a0s actuais 
diplomados pelos Institutos Industriais possibilidades de acesso à categoria de 
engenheiros ; 

realização nas universidades, nomeadamente nas escolas superiores de engenharia, 
de trabalhos de investigação, nos quais participem professores, assistentes e alunos ; 
funcionamento de cursos pós-universitários para os engenheiros, os quais, per- 
mitam manter actualizados os conhecimentos destes, tanto no que respeita às 
matérias aplicadas, como às ciências de base e na efectivação dos quais participem 
activamente as escolas e entidades ligadas à engenharia; 

activa participação da indústria na resolução dos problemas do ensino, inclusivamente 
no que se refere aos meios materiais. 


II— É indispensável a criação de organismo que defina, em função dos parâmetros econó- 
micos e sociais do País, qual a qualificação e a quantidade de mão-de-obra, em todos 
os escalões da actividades profissional, com vista a: 


a) 


h) 


c) 


, 


instituição, à escala nacional, de uma política única relativamente ao desenvolvi- 
mento economico, ao ensino e à investigação; 

estabelecimento dos escalões de nível do ensino para as diferentes actividades 
profissionais ; 

alargamento do âmbito do ensino de modo a cobrir eficazmente, de acordo com 
as necessidades do País, um maior número de ramos especializados, nomeadamente 
nos cursos de engenharia. 


WI—E fundamental que o acesso aos níveis superiores do ensino seja facultado ao maior 
número possivel de indivíduos, para o que se torna necessário: 


a) 


b) 
c) 


adopção de medidas de carácter económico — nomeadamente a concessão de bolsas 
e subsídios — que tornem possível o prosseguimento dos estudos a todos os portu- 
gueses que, em qualquer ponto do País, para tal se revelem aptos; 

criação de um sistema eficiente e generalizado de orientação profissional que 
permita fazer uma selecção adequada; 

instituição de múltiplas modalidades de ingresso nos diversos escalões do ensino, 
de modo a estimular em todos os portugueses o gosto pelo estudo e pelo progresso 
intelectual e profissional. 


IV— Considera-se primordial, dentro do problema geral da rentabilidade do ensino, a eficiência 
dos meios humanos que para ele concorrem, para o que se torna necessário: 
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a) 
b) 


constituição de quadros docentes suficientemente amplos com elevados níveis cien- 
tífico, pedagógico e cultural; 

estudo das melhores condições de recrutamento e de permanente actualização do 
corpo docente ; 


c) integração dos professores e dos alunos na política dos objectivos da Escola, 
considerando-se uns e outros como profisssionais nas respectivas actividades; 

d) existência, em cada escola, como elemento de interesse fundamental no ensino, 
de uma associação de estudantes que, devendo actuar em íntima colaboração com 
o corpo docente, contribua para o equacionamento dos problemas pedagógicos 
e para uma mais perfeita e completa formação dos futuros engenheiros ; 

e) apoio às demais organizações circum e para-escolares que, de algum modo, concorram 
para os fins enunciados na alínea anterior; 

f) criação de condições que conduzam à participação dos engenheiros nas activida- 
des das escolas. 


V— Torna-se imperioso encarar e com realismo resolver o problema das verbas necessárias 
para apetrechar as escolas de engenharia, considerando: 


a) a necessidade do equipamento laboratorial de custo elevado, para ministrar um 
ensino actualizado que acompanhe os rápidos pregressos da Ciência e da Técnica; 

b) a elevada preparação científica e técnica exigida para a utilicação desse equipamento; 

c) a necessidade de um quadro de pessoal auxiliar (preparadores, mecânicos, electri- 
cistas, etc.). 


VI — É necessário renovar a estrutura de funcionamento das escolas de engenharia, mediante : 


a) modificação da legislação vigente de modo a conferir-lhes mais autonomia adminis- 
trativa e pedagógica, dentro de uma orgânica geral a respeitar, em especial através 
da criação de órgãos consultivos nos quais estejam representadas as entidades 
públicas e privadas interessadas no ensino da engenharia; 

b) orientação do ensino com um carácter essencialmente formativo ; 

c) planeamento e coordenação dos programas e sua actualização permanente, em 
especial no que se refere à interdependência das cadeiras de preparação científica 
e da especialidade: 

d) estabelecimento de um quadro equilibrado de professores repartido desde os que 
se têm de dedicar às escolas em regime de full time, até àqueles que apenas vão 
ensinar matérias da sua especialidade, com a correspondente diversidade de 
remunerações ; 

e) apreciação do aproveitamento dos estudantes, com base no trabalho comum de 
professores e alunos. 


VII — Considera-se da maior importância a continuação dos trabalhos iniciados por este 
Congresso, através da criação de uma comissão permanente com a função de manter 
o debate sobre os problemas relativos ao ensino de engenharia e coordenar os estudos 
que, neste domínio, venham a ser empreendidos. Esta incumbência é, provisóriamente, 
atribuída à Secção de Antigos Alunos da Associação dos Estudantes do Intituto Supe- 
rior Técnico, até à criação de um organismo que englobe representantes do ensino, 
dos diplomados e das entidades ligadas à engenharia. 


Lisboa, 26 de Julho de 1963 


A COMISSÃO REDACTORA 
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RELATÓRIO DA COMISSÃO ORGANIZADORA 


Ao chegar, práticamente, ao fim da tarefa que lhe tinha sido atribuída, a Comissão Organiza- 
dora julga conveniente tornar público um pequeno resumo das suas actividades. 

A oportunidade de um congresso do ensino superior de engenharia foi verificada pela Comissão 
Organizadora das Comemorações do Cinquentenário do Instituto Superior Técnico que se propôs 
levar a efeito tal realização. Para isso escolheu uma comissão de três dos seus membros (José Maria 
Mercier Marques, Paulo de Barros e Décio Thadeu) para elabórarem o respectivo programa e, uma 
vez aprovado este, nomeou a Comissão Organizadora e a Comissão Executiva encarregadas de lhe 
dar realidade. 

A partir desse momento, toda a responsabilidade dos trabalhos ficou entregue àquelas Comis- 
sões, cuja constituição é conhecida. 

Deve, no entanto, esclarecer-se que a Comissão Organizadora nunca ficou completa, dada a 
impossibilidade que se verificou em nomear um delegado dos alunos da Faculdade de Engenharia 
da Universidade do Porto, mau grado as tentativas reiteradas que se fizeram, quer por intermédio do 
aluno José Fernando Pinharanda, delegado da AEIST, quer por intermédio dos delegados do corpo 
docente daquela Faculdade, as quais ficaram sem resposta. 

Só após a realização do Congresso, foi possível agregar à Comissão Redactora das Conclusões 
Finais um aluno da Faculdade de Engenharia, membro do Congresso, — Luis Vaz Pinto. 

Ainda durante o período de organização, verificou-se a substituição do Prof. Manuel Corrêa 
de Barros pelo Prof. Augusto Farinas de Almeida, em consequência de o primeiro ter sido nomeado 
Reitor da Universidade do Porto. 

A realização do Congresso, integrado no programa das Comemorações do Cinquentenário 
do Instituto Superior Técnico revelou-se, desde muito cedo, impraticável. 

Assim, obtido o patrocínio do Ministério da Educação Nacional, dado pelo respectivo Minis- 
tro, Senhor Prof. Francisco de Paula Leite Pinto, e renovado, posteriormente, pelo Senhor Prof. Ma- 
nuel Lopes de Almeida, procurou a Comissão Organizadora assegurar os fundos indispensáveis 
a uma organização, que embora, modesta, garantisse o mínimo de decência 

As diligências efectuadas nesse sentido, junto da Fundação Calouste Gulbenkian e das empresas 
industriais, foram laboriosas e longas ; paralelamente, surgiram circunstâncias adversas na vida nacio- 
nal que obrigaram a sucessivos adiamentos, só se podendo verificar a realização do Congresso na 
última semana de Novembro de 1962. 

Posta, por várias vezes, essa hipótese da não realização do Congresso, devido às circunstâncias 
adversas e imprevistas surgidas, foram sempre, a Comissão Organizadora e a Comissão Executiva, 
de parecer unânime que, não obstante todos os adiamentos tornados inevitáveis, a realização 
do Congresso se devia manter. 

Os adiamentos sucessivos se, por um lado, justificaram até certo ponto alguma descrença 
manifestada na organização do Congresso, permitiram por outro lado, obter a colaboração efectiva 
da Faculdade de Engenharia da Universidade do Porto que, desde início, a Comissão Organizadora 
das Comemorações do Cinquentenário do Instituto Superior Técnico julgara conveniente. 

Chegados a bom termo os trabalhos do Congresso e publicados os relatos das sessões 
de trabalho, as conclusões das secções e as conclusões finais, a Comissão Organizadora considera 
terminada a sua tarefa e tem a convicção de que a persistência posta na execução do Congresso 
do Ensino de Engenharia foi compensada pelos resultados já obtidos. 

Estes não foram excepcionais; no entanto, provaram o grande interesse de os assuntos do ensino 
serem largamente debatidos e, se as conclusões não têm força executiva, não será possível 
desconhecer, no futuro, as directrizes mais importantes apontadas. 

A Comissão Organizadora pensa que a efectivação do Congresso constitui um serviço impor- 
tante prestado, quer ao País, quer às escolas de Engenharia. 

Muito do mérito da realização do Congresso do Ensino de Engenharia não pertence à Comissão 
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Organizadora e à Comissão Executiva, embora estas tenham procurado polarizar em si todas 
as contribuições e boas vontades. 

Em primeiro lugar, deve destacar-se o importante subsídio concedido pela Fundação Calouste 
Gulbenkian que permitiu dar início público aos trabalhos de organização e sem o qual, é quase 
certo, o Congresso não se teria efectuado. 

Infelizmente tem de sé concluir que o Congreso não encontrou grande eco entre a Indústria 
Nacional. Com efeito as diligências realizadas para obter algum apoio material por parte daquela 
conduziram, na maioria dos casos, a respostas negativas (1). 

O número dos Congressistas, cerca de seiscentos e cinquenta, considera-se muito satisfatório, 
atendendo a que era a primeira vez que se realizava uma iniciativa deste género no País, facto 
que justificava alguma descrença e, sobretudo, expectativa. 

Não foi possível garantir um serviço permanente de secretaria (os recursos de que se dispunha 
não o permitiam); todavia, conseguiu-se assegurar a eficiência indispensável, quer na organização 
geral do Congresso, quer na realização das sessões, que se ficou a dever à dedicação e à competência 
do pessoal da Contabilidade e da Secretaria, bem como do Chefe do Pessoal Menor e Contínuos 
do Instituto Superior Técnico. 

O mérito do Congresso deve-se, em primeiro lugar, aos Congressistas que apresentaram comu- 
nicações, tanto mais de agradecer quanto pouco numerosos, e sem cuja contribuição o Congresso 
não teria realização possível. 

Ao trabalho dos Relatores (realizado, aliás, em precárias circunstâncias devido ao tempo 
extremamente curto que a Comissão Organizadora lhes pôde conceder) ficaram a dever-se as sínteses 
que serviram de base ao funcionamento do Congresso e a redacção das conclusões das Secções. 

Os Presidentes e os Secretários das mesas asseguraram a efectivação das sessões com elevada 
correcção e competência. 

Grupo numeroso de Congressistas participou activamente nos debates, permitindo dar-lhes 
maior generalidade e manter vivo o interesse pelos trabalhos, contribuindo decisivamente para lhes 
impor nível elevado. 

Finalmente, a Comissão Organização contou com a colaboração valiosa dos Congressistas 
Eng. Manuel Rocha, Ag.-Técnico de Eng.2 Henrique Ramos Antunes, Aluno da F. E. Luís Vaz Pinto 
e Aluno do 1.S.T. Júlio Freire de Andrade para a redacção das conclusões finais. 

A realização do Congresso do Ensino de Engenharia foi tarefa muito espinhosa, sob todos 
os aspectos, que só foi possível levar a bom termo mercê da valiosa colaboração recebida e que 
se registou com agrado. Assim, a Comissão Organizadora quer manifestar a todos sincero reconhe- 
cimento e tem o dever de lhes endossar a maior parte do sucesso alcançado pelo Congresso. 


Lisboa, 27 de Julho de 1963 


A COMISSÃO ORGANIZADORA 


(1) N. R. — Os nomes das dez empresas que contribuiram para a realização do Congresso constam das primei- 
ras páginas da publicação editada pela Comissão Organizadora. 


TEONICA 
5 


NOTAS INFORMATIVAS 


C. D. U. 624.311.5/94 


Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 


— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas (R, N. C.) 


Nora: As produções e os consumos das empresas do R, N. C. representam 
cerca de 94 º/, dos totais do Pais, 


I — Breve nota mensal 

No aspecto hidrológico, o mês de Julho apresen- 
tou-se, no conjunto, muito abaixo da média. 
[1 — Elementos gerais (GWh) 

a) Mensais 


Variação 
1962 | 1963 o | 
Ú 


SS i.e | | e | 


Produção hidráulica (Ph). 804,7 | MUT,S| + 1 
Produção térmica (P+)..... 2,1, 18,4 — 
Produção total (PT). ..... 806,8 | 826,2 | + 6 
Exportações (Ex) .. ...... 0,0 0,0 0 
Importações (1 ..o sa sys a 0,0 0,0 0 


Produção para con- 

sumos não perman. (Pcnp) 
Produção para con- 

sumos permanentes (Pcp) (!) 
Coeficiente de hidraulicidade 


693) 615) — 1 
237,0 | 264,7, +411,5(2) 
0,86 | 0,68)  — 
b) Acumulados desde 1 de Janeiro 


Variação 
1962 | 1963 0 % 


Produção hidráulica (Ph)... 2129,8/2357,7|+ M 


Produção térmica (P+)...... 82) 37,8]+355 
Produção total (PT)....... 2138,0/2395,0| + 12 
Exportações (Ex) ......... 0,0! 32,5 — 
Importações (1)... ........ 0,0! 01 — 


Produção para con- | 

sumo não perman, (Pcap) | 156,3 | 499,9 | + 10 
Produção para con- " 

sumos permanentes (Pcp) (!) [681,1 1892,4 | +-10,7 (2) 
Coeficiente de hidraulicidade | 1,13] 1,44 — 


NOTAS : 
(1) A produção para consumos permanentes (Pcp) é determinada 
pela seguinte expressão: Pcp = Pr — Penp Ko ; A DS 


(2) O aumento percentual de produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, é 
respectivamente de 10,6 º|, e 10,3 º/,» 


III — Diagramas de carga dos dias característicos 


4,* teira: 


1962 1965 

Produção hidráulica (Ph) MWh | 10430 10845 
Produção térmica ( Pe) MWh 0 0 
Produção total (Pr) MWh | 10450 10845 
Trocas com [ Export. (Ex) MWh Ú Ó 
Espanha O Eno (1) MWh U Õ 
TOTAL Pr -+(I-Ex) MWh| 10430 10845 


Prod. para cons. perm.  (Pcp) MWh| 8452 9342 
Prod. para cons, não perm, (Penp) uwh! 1978 1503 


Potência máx. MW | Ra 608 — 
ads ta ; Potência min, MW 10,09 284 
SE id Utiliz. da ponta horas 18,1 17,8 
ER Factor de carga Oo 0,74 
“E [Potência máx, MW) 49 540 
E q p Potência min. MW 189 215 
ós cp | Utiliz, da ponta horas 17,3 3 
Factor de carga 0,72 0,72 
TRONICA 
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IV — Energia armazenada nas principais albufeiras 


| | No fim do mês 
Albufeiras : 

GWh 9a (4) 
Paradela. . . . cc... «| 188,9 | 62,4 
Venda Nova +. «. cv... 79,6 62,2 
Salamonde . . «cce 22,8 62,6 
DRC emo e sas sw» 21,4 64,7 
Cabrili .. 20... «| 2869 81,6 
Castelo do Bode. . ..,..| 1038 63,4 
Guilhofrei «mw pon mw ou TA 89,2 
Lagoa Comprida . . +. ... 27,1 (2) SL,O 
Santã Luziã . ccc move õ1,4 83,4 
PrACANA pcs russo 1,9 61,2 
PÓVOA náo cs sz 12,9 (3) 98,1 
Total. . 199,6 12,8 

Notas: 


(1) Coeficiente de enchimento em energia das albufeiras, 


(?) Inclui 3,8 GWh armazenados em Vale do Rossim no 
fim do mês, 


(º) Inclui 2,7 GWh armazenados no açude do Poio no 
fim do mês, 
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Tratamento matemático dos movimentos vibratórios 
de características lineares 


PELO ENG.º MEC. (1.8. T.) JORGE NEVES DA SILVA 


(Assistente no 1.5. T. — Bolseiro do 1. A. €. 


o. INTRODUÇÃO 


Neste trabalho pretendemos apresentar as soluções dos modelos matemáticos construídos a 
partir dos modelos físicos idealizados de sistemas vibratórios reais. É corrente na solução destes 
modelos apresentar soluções e verificá-las, introduzindo coeficientes que se definirão na substituição. 
Como veremos, é possível deduzir as soluções integrando com aproximações sucessivas e determi- 
nar imediatamente os coeficientes com toda a generalidade. 


1. SISTEMAS COM UM GRAU DE LIBERDADE 
1.1. VIBRAÇÃO LIVRE AMORTECIDA OU NÃO 


1.1.1. Modelo matemático 


Consideram-se aqui os sistemas constituídos por uma massa sujeita a um movimento simples 
linear, actuada por uma ou mais molas de qualquer tipo, mas de característica linear, dissipando ou 
não energia em atritos de carácter viscoso (*), e livre de quaisquer acções exteriores. Poderemos 
substituir todas as molas por uma equivalente de característica k, e os atritos por um amortecedor 
linear equivalente e de característica R. O modelo físico resultante será, fig. 1.1a ou fig. 1.1b. 


Á 


x(t) k R 
m 
x (t) 
Fig. 1.1a Fig. 1.1b 

o modelo matemático correspondente (**) será 

dº x dx 

m + R— + ki==0 

dt dt 
ou 

x: + Ed x + EA x==0 

m m 


(*) Entende-se por viscoso um atrito em que a força resistente varia linearmente com a velocidade. 
(**) Ver, por exemplo, À. Portela — Máquinas Alternativas I, 2.º Livro. Ed. A. E. I.S.T. 
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R k 
Usando a notação corrente, 2 = E e Wo = a 4 respectivamente factor de amortecimento 


e frequência própria do sistema não amortecido, será 


x+2)Ix+to0)x=0 (1) 
No caso particular do movimento não amortecido, R = O e 
x 0,3 x = O (2) 
1.1.2. Integração da equação diferencial x + 20 x + wlx=0 
Fazendo 
xXx == %4 
XxX = = x 
x=u=-—20x — 0,2 X 


== — 20x3— 0, 4 
a equação matricial 


ou, usando uma notação breve, 
DX = aX (3) 


representa as expressões acima consideradas. 
A matriz quadrada « tal como se apresenta é de difícil tratamento nas operações que se vão 
seguir. Convém efectuar uma mudança de coordenadas orientando os eixos segundo os vectores 


próprios de a: 
TºDX=TAasTiX= ATX (4) 


Determinação da matriz diagonal A : 


det [x— AI] =0 (5) 


det ua Ê . = 0 
— to? sm DO cmi À 


»[29+1] + o? =o0 


ou 


equação cujas soluções são 
odio +iVo—-?=—)+-io 


a=—)—iVoi-d =—j—ia 


(6) 


usando a notação usual, em que w representa a frequência própria do sistema amortecido. 


MH O —)+io 0 
A = sum 
o ps! | O 4 (7) 


Obtém-se a matriz de transformação do seguinte modo : 


Teremos então 


T-* aT= A 
&T= TA 
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Tendo em atenção os vectores próprios da matriz «: 


[a-M 1) | (é |=0 ou ta =D tu 
tai 
[ami] | 2 ]=0 ou tag == do tua 
taz 
será 
tu tis a MTE SE IA. | t 1 
T= = tu t = tuts 8 
o fed] - a), 1 | + io Media (8) 
e 


Ties La erasal cr alo mf 
tibs (0a — 9) L—dl 1 2 o tati —io 1 
Voltando à equação (4) e fazendo Y = T-!X, 
DY AY (9) 


A integração desta expressão poderá realizar-se por aproximações sucessivas: 
t t t t é 
Y=Y+ | NAYd=Y+ A [Yo +f A|Yo+ A (Yo + n(::-) de) de | de [at= 
to to to to io as 


= [14 ot) 4 nt tod po ir 4 ..ly, 


Esta série corresponde ao desenvolvimento da função exponencial. Portanto 


A (t—to) 0 
[v] = et) [ygj= [Ye] (10) 
As (t ui to) 
0 e 
uma vez que [A] é uma matriz de elementos constantes e é diagonal. 
Efectuando agora a transformação inversa 
Ty TIA AX o X 
e 4 (t— to) 0 
X = T ET NY 
0 é da (t — to) 
[A] (t—to) qe 1 *. & M(t — to) o à —1 
aa AR aÃ Rd 
da— hM [LM dy o e — M 1 
a ES to) MU St — to) — et — to) 4: Et — t.) cú 
Àg — hM MM ha im to) cs M Às aut to) na | tt sir to) fo Ja tt Má to) ' 
Como T Y, == X,, será, substituindo À e 72 
E A ; ; e. 
X1 cos[w(t—to)| + = sen[w(t — to] e sen[w(t — to)] X04 
=. S(t — to) i 
x2 — - sen [w(t — to)| cos[w(t—to)] — FA sen[w(t — to)] x02 
(12) 
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Finalmente, tendo em consideração as substituições efectuadas no início e fazendo x,, = Xo 


O Xw =X, vitã 


— 3(t— to) 3 + x 
es E | cos[m(t—t,) + eo 


sen [m(t — to)] | 
. (13) 
y — S(t — f sa Xa 0 Xo R 
A to) |— Eee R sen[o(t—t,)| + x, cosfo(é—t5)] | 

A primeira equação de (13), integral geral da equação homogénea (1), representa o movimento 
da massa em função do tempo. A segunda, que aliás se poderia obter directamente da anterior por 
derivação, representa a velocidade da referida massa. 

Confirmam-se as definições de » e à referidas atrás pela observação da expressão (13), pois 
em *(t—to) é um factor de amortecimento do movimento enquanto que w exprime uma velocidade 
angular. 


Se não existirem amortecimentos, ) = 0 e 


= x, cos [m, (t —to)] + Fe sen [ou (t—t )] (14) 


X = — 0 X, sen [mo (t—to)] + Xo cos [5 (t— to)] 


Estas expressões poder-se-iam determinar directamente a partir da equação (2) por este ou 
outro processo (*). 


1.2. VIBRAÇÃO FORÇADA AMORTECIDA OU NÃO 
1.2.1. Modelo matemático 


Neste estudo consideram-se só as perturbações de carácter sinusoidal introduzidas nos sistemas 
cujo modelo físico se apresentou em 1.1., o que equivale a considerar qualquer perturbação periódica, 
uma vez que sendo esta decomponível numa série de Fourier e dadas as características de linearidade 


(*) Um processo de integração da equação x + wo? x == O que apresenta amplas perspectivas no estudo de 
movimentos não lineares, será estoutro : 


xxtu?xx—o 


x + (dg? x* == x; + wo? %o* = 2 U 


ou 


exprimindo U, a energia total do sistema em qualquer instante, valor constante uma vez que não há dissipação 


dx —— V2 U, us (9,2 x 
dt 
"K dx dd 
| E E dt 
V2 Us, — tw? x 
x, to 
ou 
E = 
1 V2 U, Som? 1 tio x 
Tp ———açem ; "| arc cos /2U x = [— E 
tr LO BR tos 
ve SM Val, da 


ou, finalmente 


/ 


V xo? + w2x? 
Xo 


Va? Ag 


o E (t — t,) + arc cos 


Vo 


== Xo COS [0 (t — to)) + — sen [o it — to) ] 
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exigidas aos modelos elásticos, o movimento resultante da perturbação complexa será um somatório 
de movimentos resultantes das perturbações sinusoidais de frequências diferentes. 

A perturbação F (t) pode actuar na massa como força (fig. 1.2a) ou nas ligações ao exterior 
como deslocamento, (fig. 1.2b) mas sempre na direcção do movimento. 


x (t) 
m 


F(t) 
Fig. 1.2 a Fig. 1.2 b 


Teremos para modelo matemático (*) no primeiro caso : 


mx + Rx +kx=EF(t)= z | Ain cos [nB (t — to)] + Agn sen [nB (t —to)] | (15) 


No segundo caso, consoante o movimento perturbador actua parcialmente ou totalmente as 
ligações, assim teremos uma equação (15) em que os coeficientes A; representam amplitudes de 
movimentos em vez de amplitudes de forças, a menos um factor de dimensionamento. 

Tendo em consideração o que se disse no início desta alínea bastará esquematizar o problema 
para uma perturbação sinusoidal pura, por exemplo, a de frequência fundamental. 


x+23x + 02x= É sen [B (t—t))] (16) 
m 


ou, para um sistema não amortecido )=0 e 


cito q gu “a ses [É tes) (16') 


O integral geral para esta equação diferencial será 
x = xg + xp 


em que xy significa o integral geral da equação homogénea ou seja em regime livre, e xp O inte- 
gral particular para o sistema forçado. É este último valor que falta determinar. 


1.2.2. Integração da equação diferencial 


ds cal E 
x+20x+wx=—sen[B(t—t)] 


Tal como em 1.1.2. virá 


O 
Fo! 


DX = ax + “a a Bo | (17) 
n E = aX + Xp 
(*) Ver, por exemplo, A. Portela — Máquinas Alternativas | — Ed. A. E. 1.5. T. 
TRONICA 


A equação (4) ficará 


[DX TtaTTAxX + T+HXr (18) 
sendo agora 
Yp== T+ Xp me To | É | me sen [B (t—to)] = [A] sen [B (t—tç)] (18) 
E a equação (9) fica 
DY = AY + Ye (19) 


A integração, aplicando o mesmo método, será 
v=y, +['vedt+ [nva=y, + [Pratt Alvo | ve att [PA DrA+ | Ye de + 
to to to to to to to 


+ ne [ve a + [ A (...) dt) de at di 
to to ” 


A parte correspondente ao integral particular virá, depois de substituir (18”) 


Yo= E Laje [e [ [4] sen [B(t—t)] de = 


= [— 1) cos 8 (t=t,)] — “Al sen [B(t—t9] 4 IG-8 


| [AP O, LAS 1 A] 
DITAS. co BEMA qua 30) SOS 
Li 
e, uma vez que 
[AJt = [Aiitl 
por ser uma matriz diagonal, 


[A] [0] (tt) 


ze | — I cos [B (t—t5)] — + sen [B(t-t)] +e 


a 


TORA DEN (8 o 
ce relé) 


Este desenvolvimento só será válido se as séries dos hii forem convergentes. É condição 
suficiente de convergência absoluta que 


mod 4 << B 

mod )2) << B 
ou seja 

0! + ut<B 
ou, uma vez que wi =m9?-—-0d, wy<B 
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Nestas condições 


Yp = |—Icos [B(t—t)] — EA sen [B(t—t,)] + A —t) 


LA] = 
: (0) id (2) | 


— cos [B (t—t)]— sen IB (t—t)] + e (t— to) 


—— 


a O 
0 E 
da Aa (t — to) 
— cos [B tolo sen |[B(t—t)l+e 
go 
Bº + A? B [A] = 
Es Bº + )9? 1 
Sa wÊ 
=| 8 2 | ta1=tt4 e 
em que 
À da bom 
pá = ET) |—cos [B (t — t,)] =. sen [B (t—to)] + e” t é 
jm ga |—cos [B (t—to)] — x sen [B (t—to)] + di 
Bº + Ao? B 


Efectuando agora a transformação inversa 


= TL=TyT ITA (22) 
Da equação (18º) 
o ros = a = 
O 
Fos 
E + — 4 es T— 
TA=TI E Fo 


Tendo em consideração (8) e (8”) 


da —1 | CMq—A ” ? 
Tot t=[ oil pa fa |=+ 2 eu da pa 1 + pa | em 


20 LM JILO qall—-m 1 20 Lhda(o—ç) —Mp + J202 
Finalmente 
” : ha 01 — My qa — 1 + 02 » i Eu —p1 + fa | 
RT Aa 01 ZE 
a” da da (04 — 09) —he + da = em — ut + 2202 
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e substituindo valores para /M, /3, o e &3 e tendo em atenção que 


Xp a Xi p | ai Xp 
- X2p - Xp 
| Eos 
fd PP | 20) Bcos [B (t—t,)] + (Bº--,?) sen [B (t—t,)] — 
| 2 ras 
-— ad E (BP+2)—m,2)sen [w(t—t,)] + 238 cos [o (t—49]]] 
Fy 
Xp = — CP PRE: | — 205Bº sen [B(t—-t)]+(B'—-w?) Bcos [B(t—t,)] + 


a 


dig [E (B? + 05º) sen [o (t—t,)]- B (Bi 0.) cos [9 (t—,)] | (24) 


Estas expressões representam o deslocamento e a velocidade do sistema com toda a generali- 


dade desde a posição de repouso. 
No caso de não haver amortecimento, )=0 e 


| Fo 
B 
Xp == — RR 1 sen [B (t—-to)] — o sen [mo (t—to)] | 
| pão as) 
m 
[Xp = BiDos | cos [B (t—to)|] — cos [o (t—t9] | 


expressões que se poderiam deduzir directamente a partir da expressão (16). 
A solução completa para a equação (16) será então, adicionando as expressões do integral 
geral da equação homogénea, e do integral particular agora determinado 


e B (Bº — 6?) For 


| =St= 1] | 25 B ot 1 
ii “o + snes e 1 Baal io 
dai B+28 o) ê 1 
o“ o ; / 
Tl * (Bi odprapra | bo tt—to)] | já 

Fo: 
m , 

| -—— (BD o? E TF | 20)Bcos [B(t—to)] + (B'—-«ç?) sen [B (t—to)] | 

1 (26) 


e at (t—to) ] 
x= e Xo + ——— o ns 
| = pano [e be fab | 
B3 S » Fo 
et xe , Co Pr - 
[od Ceasa |ºo to tod] [— 
Fo 
ai RES == | (8! — 00%) B cos [B (t—t)]= 22 B' sen [B (t—t9]| 
| (Bº — wo”)? + 4 Bº 3º | 
TEONICA 
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“Perbunan C - evita a chama e resiste ao calor 


Os artigos feitos da borracha de cloroprene Perbunan C, distinguem-se por uma combinação favorável de 
propriedades tecnicamente valiosas. Evitam a chama e resistem ao calor, às intempéries e ao ozono, aos óleos e 
aos agentes quimicos e têm boas propriedades mecânicas. 

O Perbunan C é especialmente apropriado para a fabricação de bainhas de cabos eléctricos e isolamentos para 
corrente de baixa tensão, artigos técnicos moldados, para a construção de automóveis, de vagões e de máquinas, 
para perfis de vedação, mangueiras, correias de transporte, rolos, para a impregnação e o revestimento de tecidos, 
para solas e para colas. 


PERBUNAN C o Cloroprene da BAYER 


Representante em Portugal: 

S.A.R.L. QUIMICOR 

Lisboa: Rua Sociedade Farmacêutica, 3 
Tel. 555496 e 555893 


2 Porto: Rua de Santos Pousada, 441 N » 
çã Tel. 54141 P SS. o TlialoeTla. 1863 q» 1963 
DA 
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Sobre 2 ou sobre 4 rodas 
Peças de plástico Cellidor 


O “Cellidor é um excelente material para a 
fabricação de acessórios para automóveis e para 
biciclos. Não tem apenas boa resistência mecânica, 
também as suas superfícies conservam e 
aumentam o seu alto brilho, pelo uso, e tem 

ainda um toque agradável. 


Presta-se também para envolver peças 
metálicas, pois não abre fendas, que dariam 
lugar à corrosão. 


Este material amortece o som. 


Não atrai as poeiras porque não é susceptível de 
carregar-se de electricidade estática. 


Há muito que o Cellidor se introduziu com 
vantagem ne fabricação de acessórios para 
automóveis, como volantes, revestimentos das 
colunas de direcção, botões de comando 

nos «tabliers», revestimentos para 
conta-quilómetros, farolins trazeiros, frisos 
interiores, assim como protectores de 
correntes e chapas protectoras para biciclos. 


Quando precisar de tipos especiais para os seus 
artigos, consulte as nossas representações 


S.A.R.L. QUIMICOR, Representante em Portugal 
Lisboa, Rua Sociedade Farmacêutica, 3, Tel. 42194 
Porto, Rua de Santos Pousada, 441, Tel. 54141 


E 


B * 


c> el Ê F al €> me ses (bs jones 
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Ficam assim definidas por via analítica as expressões completas do deslocamento e da velo- 
cidade do sistema perturbado com condições iniciais quaisquer, englobando as fases transitória 
e permanente do movimento. 

Ou, para o sistema não amortecido 


Í Es E 
RE da a 
x =xo cos [wo (t—to)] + = + a (Bi 0,3) sen [mo (t— to)] — Bios Sen [B (t — to)] 
8 j 
- B For = B Fo 
m m 
| Xo = | BiDoê + xo | cos [m(t+t)] — w,x,sen [mn (t—t;)] — BiDod cos [B(t—to)] 


(27) 


Note-se que estas expressões foram deduzidas para o caso de ser B >»,, pois só nestas con- 
dições provámos a condição suficiente de absoluta convergência da série de integrais, porém esta 
circunstância não invalida a aplicação dos resultados ao caso de ser B<.», porque as expressões 
(24) verificam a equação diferencial sem quaisquer restrições paramétricas. A solução não é incom- 
pleta porque o integral da equação (16) é combinação linear de duas soluções particulares linear- 
mente independentes. 


1.3 REPRESENTAÇÕES GRÁFICAS 
1.3.1 Vibração livre 
1.3.1.1 Sem amortecimento 


Nos sistemas não amortecidos, tal como as expressões (14) mostram, a representação do 
movimento num diagrama deslocação ou velocidade-tempo é uma sinusoidal pura. Arranjando as 
referidas fórmulas 


2d) 
Xn 

x == Vas di cos [mo (t—t)— +] 
o 


(14º) 
x = Vas 0,2 + x, sen [mo (t—t,) —:] 
em que 
tg e= 
5 OR X, 
A sua representação gráfica será (fig. 1.3). 
1.3.1.2 Com amortecimento 
Nos sistemas amortecidos, e arranjando as fórmulas (13): 
a 2 a h) (t = t,) 
Xn9 + x, 
E Vas + eia) - cos (w(t—to)— 4] 
= ani (13º) 
. Ff (1,2 X ix 3 —S(t—ts) 
Xx=— V Xo? + (eu-xerde) e sen [m(t—to) —=a] 
em que 
ti di Es Xo 0 + Xe 
hi w) Xo 
e ; 
XxX, 
tg é = — Rs a 
WD? Xo + 0 Xo 
TRONICA 


Representam curvas periódicas de carácter sinusoidal de amplitudes decrescendo logaritmica- 
mente (Fig. 1.4). 


Fig. 1.3 Fig. 1.4 


Verifica-se que, teóricamente, o sistema estaciona ao fim de um tempo infinito. Na realidade 
havendo atritos de carácter não viscoso, o sistema pára num tempo finito. 


1.3.2 Vibração forçada. 
1.3.2.1 Sem amortecimento. 


As expressões de (27) poderão escrever-se 


Eu 
- note É om fig email Ric = 
— Xo RE cos [o, (t—t;)—: à (B-+r0,) sen [wo (t—t, 
Fa 
m (B—m,) (t—t,) (B-+o,) (t—t,) 27º 
nc B'— o? sem RR cos E ( ) 
nu PP 
CREA A m — 9) (t — 5) (t—t, 
XX == — Vs we: x, senjo,(t—-t)—:] + ei sen a — a — RR 
a É. 
em que tg * = DX; 


Pela expressão da velocidade verifica-se que esta é dada pela sobreposição de um termo sinu- 


à Ê (1). " 
soidal de frequência —” e de um outro representando um «batimento» de wo com B. 
2% 


Se no momento de aplicação da perturbação for x, == xo = 0, a curva da velocidade será (Fig. 1.5). 
Nestas condições a curva do deslocamento não é tão pura, pois é uma soma duma curva deste 


º * . A . 6) - . A - + . 
tipo com uma sinusoidal de frequência e (fig. 1.6), sendo a importância desta última curva tanto 


menor quanto mais se aproxima B de wo. 
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Fig. 1.5 Fig. 1.6 


No caso particular de ser B = wo, portanto quando a frequência imposta ao sistema entra em 
ressonância com a sua frequência própria 


des B-o)(t>t) — (B+o)(t-to) 
= ce— (o) — 
lim >. sen LE NES o 
B* — mp? 2 2 
B — og 
Fou  B-m “ 
— sen | ———— (t—to) 
m e - Jo. (B-+-mo) (t—to) 
lim Pr] Rs (t — to) cos ra 
B -— Wo ' e (t —to) 
Fou 

B + t—t 

e a = (t — to) cos E em dE e) 0) 

e as expressões (27) vêm 
Fog 
x0? m 
li é + sa cos [vo (t— to) — E] + E sen [mo (t — to)] — 
Fa 
ss Sp (t — to) cos [mo (t—to)] ni. 
2 00 : ' 
Fu 
: ————— m 
(* = — Vxo? mo? + xo? sen [vo (t—to) — E] + ES cá (t — to) sen [mo (t — to)] 


a que correspondem as curvas das fig. 1.7 e 1.8 traçadas para x, == xo = O. Estas curvas justificam 


Fig. 1.7 
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porque o sistema não amortecido não é subitamente destruído ao entrar em ressonância como nos 
faz crer o formulário correntemente adoptado. 

Há sempre um intervalo de tempo necessário para que a vibração do sistema ultrapasse os 
limites admissíveis para a sua integridade, dependente como se verifica da frequência da oscilação 


À Eu 
e do cociente —— . 
m 


1.3.2.2 Com amortecimento 


A primeira expressão de (26) poderá escrever-se 


, 


ça (Dai p x é 2 
nm V x3+ (uid) cos [»(t—t,) — 4] — 


ts / 
Fou 
.B EA - —— cos [u(t—t;) — :3] — 
tn) V(Bº cm i)0 8 )* + 4 9º Bº 
EM (26) 


EPP. - cos [B (t — t,) — 0] 


V(B'— mo?)? + 4 3º Bº 


em que 
tg 4 = x, 9+*, 
0) X» 
k Bº + 29º e, as 
gt == 
20 
Dag À 
tg 6 e E (uy 
2ºB 


4 é o chamado ângulo de fase entre a força perturbadora e o movimento. 
A expressão da velocidade será 


' — 3 (t—to) - PR E ER 
X% = e Vs? + (roxo ta) cos [o (t—to)—:4] + 
Fo; 
B q — 3 (t — to) 7 
e dee =ed >> cos [0 (t= to) — 02] — 
E V(B— wo)? + 4 3º Bº 
Fo, 
m (26") 


—E — — cos [o (t—to) —t] 
V (B? — mo?)? + 4 5º B? 


em que 
E Wo? Xo + 0 Xo 
p= 
(O Xo 
E, Õ (B2 ep wo?) 
tg Eg = — 
w (Bº — wo?) 
; 20 B 1 
Bº — wo? tg 0 
TRONIGA 


Finalmente as expressões referidas em (26) tomarão a seguinte forma: 


— & (t — to) : Re 
x=e V xo 4 (toi xe) cos [w (t— to) — :2] + 


to 


For 
DES ca ii B ed | 
Rs o (t—to)—:] — 
V (Bº— mo?) + 4 52º B! | td cos [mw ( 0) 2] 
á Fo, 
=... aa OD) De ti ia 1 BAD) e de 
“va (Bº — wo?) + 492 B? 2 2 
(26º) 
o (RT 
x = € ed TR 2 ) cos [o (t— to) —:4] + 
B o , o 
ny aa 8 — t, ' 
+ re A E pls cos [0 (t—to) — 22] — 
— to)! + 4 0 a 
B Fo 
E. q o Di io AR (iq qd) ma im 
V(Bº — 93)? + 4 52 Bº 2 2 


O 1.º termo nas duas expressões de x e x refere-se ao movimento livre. 

O 2.º termo, de amplitude constante, de frequência igual à frequência prórpria do sistema, 
poderá ser nulo em relação ao deslocamento para B==w ou muito pequeno no início do movi- 
mento mas tende para um valor estacionário quando t +» o, 

O 3.º termo representa o batimento, independente do tempo, mas é influenciável pelo 2.º termo 
à medida que o tempo decorre acabando por se anularem as oscilações de frequência » e ficando o 
sistema a oscilar segundo a seguinte expressão que define o movimento depois da fase transitória. 


| For 
m 
--—— cos( B (t— to) — 0] = 
V(Bº — qo?) + 492 Bº 
Fo; 
m 
— E aee cos [B(t—t,) — 94] 
1 V(B:— 0,3)! + 40º Bº 
28 
| Es (28) 
x=— B caco RE cos|B (t—t,) — 9] 
V(Bº — 0?) — to!) + 4 9? Bº 
Ex 
m 
E a sen|B(t—to) — 04] 
| VB 0) + a dB 


portanto, com a frequência imposta, com amplitude constante mesmo para B == mo, e com um des- 
fasamento em relação à força perturbadora: 


Estas são as expressões conhecidas e correntemente apresentadas nos livros de texto. 
Interessa apresentar gráficos das amplitudes do movimento depois da fase transitória, e do 
respectivo ângulo de fase em função 
dos vários parâmetros definidores do 


sistema. | 
Para tal convém definir as re- h =0 
lações f | 
Vr o, 
a | 
Wo | 
B ? 
fp e E | 
0 M | 
x X, x , 
[3 — RSS o A | / os 


o) VE e 
respectivamente índice de amorteci- Fig. 1.9 
mento, índice de ressonância e factor 
de amplificação, sendo xest a deflexão 9; 
estática da mola sob a acção de uma 
força igual à amplitude da força per- MN di 
turbadora. - 
Teremos então 
Xo 1 
DO OS e see E o 
Xest | Vd-d)+ra4avre RR u4 e 
À (29) 
2Y4€ =0,2 
al = 
a que correspondem as famílias de 2 
curvas das figs. 1.9 e 1.10. 
Fig. 1.10 


1.4 ENERGIAS POSTAS EM JOGO 
1.4.1 Sistema livre amortecido ou não 
Na expressão (1) podemos fazer 
mxx+Rx!+ kxx=0 


e primitivando entre dois instantes to e t 


Eai f à *"t 
m 2 + q Rx dt + k El=o (30) 


O primeiro e o terceiro termos representam, respectivamente, a energia cinética e potencial 


do sistema. O segundo representa a energia dissipada por atrito, e obtém-se substituindo x conhe- 
cido atrás 


E a t —23(t—to) , é . 
foriat=R e [( tdo ) sen" [w(t —t)] + xo” cos! [w(t—to)] — 


“ to tw) 
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amas xo sen [9 (t—t,)] . cos Lo (t-t)]| dt = 


2 3 Ko)? — 23 (t— to) 
Ro (io! xo 3 xo)! [- 1 mi E o e 


2 0,º 2 w 
3 —25(t—to) , 1 — 25 (t—to) 
nad senê To (t—t)] + 57 (1-. ) |+ 
A 2 92 / — 23(t—to —23(t— to) 
+ R x; Jo a A e +— é A sen [20(t—t,))] + 
0,2 0 4 0,2 
8 — 23(t— to) xn?tLix 
z a A 
+ A e sen? [w (t—t,)]-2R e io 
— 23(t—t, 5 — 25 (t—to) 
k = e "ad La fd) — PER e E sen [2o(t—to)] + 
aço E a 
tw) (1—. 2 (t |] 
4 wo? 
No fim de n períodos 
29(t—-t)=20nT 
o(t—t) =noT=27n 
e 
dd ont 
f R$ de= O [mt xt + xM[1—e | 
o 
e, dado que 0 = a 
2m 
— 23nT 
f ES im ES t+ Dasl[1-e CRE | = 


[Ee gls 


A energia dissipada ao fim de um tempo infinito, quando o sistema atinge o repouso, será 


oo 
e 1 1 . 
f Rx dt= —kx?+-—mx,? (32) 
2 2 
o 
evidentemente, a energia total do sistema. 
Se o sistema for não dissipativo 
nT 
f Rxº dt=0 
o 
o sistema não para mais. 
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1.4.2 Sistema forçado 
1.4.2.1 Não amortecido 


A expressão (16) pode escrever-se 
mxx + kxx = Fo, sen [B (t—t)| x 


e, integrando entre dois instantes, 


X 


x* 3 
2 


nã 


t á t t ; 
af [k = f. Fo sen [B(t—t,) | x dt (33) 


Substituindo agora x dado por (27), no 2.º membro da expressão, fica 


Foi 
) Fo, sen [B (t—to)] io dt= Eu | BID, xa )J cos [oo (t—to)].sen [B (t—to)] dt — 
mm A ci 


to to 


E 
.t B = 


r atdiia 4 
a | sen (oa (t—to) | sen [B (t—to)] dt — BD mê [ 2 sen [2B (t— to) | dt | — 
“to 


to 


“BB For ” E 


= Fo || = be do | > sen [(B +00) (et Jde+ | Lsen (Boo) (tt) ] dt » 


um B: —— Wo? 
to 


Ts Do Xo 


, f + cos [(B— oo) (t—to) ] dt | E cos [(B + wo) (t— to) | at | is 
” a 2 E! 


" tô to 


2 ,.8 E sen [2 B (tl dt | = 


p +é 
RE m l 1 > 7 
= E | E | + Xo “28 Fo) cos | (B + mo) (t-tgI=1] cá 
1 , 1 
net ateste B — wo = O) | == MG: | = à = (6 — to ac 
Ta (cos [( ) (t— to) | 1)] x F EB) sen | (B ) (t—to) | 
pg 
— qui e Bro 6-)] | + ps tg (e [2B(t-t)]-1 | 
2 (B + wo) : E Bº'— wo! 4B Re ea | (54) 


Quando a diferença (t — to) for tal que 
(B + wo) (t— to) = tm. 2% 
(B— wo) (t—to)=n.27= 


em que my e nm: são inteiros, a energia acumulada no sistema desaparece, uma vez que 2B (t—to) 
por consequência também é multiplo de 2x7, e a integração acima obtida anula-se nesse intervalo. 
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JAYME DA COSTA, L” 


MECÂNICA E ELECTRICIDADE EM TODAS AS APLICAÇÕES 


PORTO LISBOA LUANDA 
Praça da Batalha, 12 Rua dos Correeiros, 14 Cx. Postal 5174 
Telef. 22812 (PPC) Telef. 327035 (PPC) Telef. 4743 


q tel e " Ea “el | 


x nd eis , ! 


Meterial fabricado nas nossas oficinas, pronto para expedição 


Fabricação de aparelhagem eléctrica de alta e baixa tensão 
Equipamentos electromecânicos completos para centrais e fábricas 
Postos de transformação tipo monobloco 

Quadros blindados, capsulados e em armário 

firmaduras de iluminação de todos os tipos 

Contadores eléctricos e instrumentos de medida 

Telemedida e telecomando 

Instalações elevatórias de água 

Instalações de saneamento e esgotos 

Ferramentas de corte 

Motores Diesel industriais e maritimos 

Fogões eléctricos 

Frigoríficos 

Ventilação e condicionamento de ar 

Aparelhagem de Raios X e electromedicina TECNICA X IX 


“tiltros 


ALGUMAS DAS PRINCIPAIS APLICAÇÕES: 


CONDENSADORES ASEA 


PARA 
MELHORAMENTO 
DO | 


FACTOR 
DE 
POTÊNCIA 


PAGA POR CONSUMO 
DE ENERGIA REACTIVA... 


SE A SUA INSTALAÇÃO 
= ELÉCTRICA ESTA SO- 
E BRECARREGADA... 


fá OP rã, 
ÉS Ê 
: mA 
Wu 
W 


Bateria de condensadores tipo CLB, 32 kVAr 


«Consulte: 


ASEA ELECTRICA L: 


SEDE EM LISBOA Eua TELES FILIAL NO PORTO 
E, RODE ARTILMANA UM, tok Uo-no Credo RDO CAMPO ALEGRE, Mk 
= tecto rrr —sesonjo/9  ASENORT Telefone PRC. — 6210617 


Água e Açúcar + Óleos vegetais, 
animais € minerais e Cerveja 
Sabão liquido «Resinas naturais 


e sintéticas e Vinhos e Vinagre ” al 


Glicerina e Tintas « Vernizes FILTROS NACIONAIS PARA TODOS cai 
Sumos de frutos » Xaropes OS PROBLEMAS DE FILTRAÇÃO 


Azeite + Industria quimica 


Inteiramente fabricados em Portugal, por 
técnicos e operários portugueses, segundo 
a patente port.n. 29.917 


pf orórre 


Dezenas de instalações a funcionar no país 
atestam a sua excelente qualidade. 


Epp ue 


E mo 
E f 
- a od =» +98 
' - MN Tiara j um. ri 
To Ru q es od re a! 
! + dsH " ) . a! Ddt de! 
4 sra! ; Bios Deo 
| Ed per so! no- $; 
7 ' ! "1 tt L -+ j o 
8 
1 ARTS H 4 ? - í- Mid 


RUA CAPITÃO FILIPE DE SOUSA, 128 
CALDAS DA RAINHA 


CRU 
ea 


E 
s 


+ y 4 

findi “pr 
e é” assita! 
E tanir, 
RUA) 
Fist d 


Ora ni >> nº são inteiros, e o menor intervalo corresponde a n2==1, e portanto o sistera 
acumula e elimina energia perturbadora entre dois batimentos mantendo o balanço total constante 


e igual à energia do sistema antes de perturbado. 
Se B=»0,, t—t,=00 e o sistema só acumula energia, como se verifica em 1.3.2.1. 


1.4.2.2 Amortecido 
A expressão (16) virá agora 


mxx + Rx! + kxx=Fosen [B (t—t,)]x 
e integrando 


x? “2 x? t “t É at ' 
mé] + [k | +) Rxº dt=— | Fo sen [B (t—t,)] x dt 
2 to 2 to É 

to o 


ou 


U—U= = [m (x2 — x?) + k (xº — e | = f. [Fu sen IB (t— t)])—R | x dt (35) 
to 


Substituindo x obtido em (26) e fazendo na fórmula resultante 


23 B? Fu 
m 
(BE — T+ apr OA 
B (Bº — 02)? Pt 
ee ua RT 
RR TT) =, 
(Bº — 9?)? + 452 B (36) 
tg? xo + à xo = E 
tw) 
B? (B? + 642) Fo 
(8) m us G 


(B'—- o?) + 452 B2. 
obteremos 


U-w=[. (Fu RA) sen [B (t — to) ] Se didi 


ia (xo + D) . cos [o (t— to)] + 


5 (t 


—— Re. 


— S(t—to) . “a 
k [e (xo + D) cos |» (t—t)]— e 


— ty) 
(E + G)sen [o (t t)] + RD cos !B (et) . 


id (E + G) sen [m (t— tJ)] — 


— Dcos B(t—-ty)]+A sen |B (t—to)|] | dt == 


S(t — to) 


e (Fu— 2 RA) (xo + D) e pm (B+o) (t—to)] + sen [(B—s). (8) | dim 


U J 


cos (B—w) (t—to)]— cos (B-+-0) (t— to)] | dt — 


— 


b) (t ei tu) 


tt 3 = 
ns di = (Fau—-2 RA) (E + G)e 


t 
e: e (Fu —-2RA) D sen (2B (t—te)] dt + 
0 
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et nt 
+ | (For A + RDº—R A? sen? [B (t—to)] dt — | R D? dt — 


to to 
at =. à EE cê 
—| RG +D-RE+GI eos [o (t—to)] dt — 
to 
ft o 3 A at R o. pas 3 
—[ rE+Ge CCT + | R GD) (E4+G) e Tt sen [20 (t—to)] dt + 
to to 
do -=—3(t-—t,) : 
+ | R (xo+ DD e [cos ((B+o)(t—to)]|+cos [(B—w) (t—to)]] dt — 


to 


— 3 (t— to) 


a R (E+G)D e [sen [(B+) (t—to)]—sen [(B— q) (t—to)]] dt = 


to 


= |[— 5 Uu—2AR) Go +D) + R (E+G)D | RG +D)D+ 


Broa 


+ 5 (248) E+O | = pot 2 sen [(B 4 1) (t—to)] + 


+ 0)P 6º 


À 7 B+o 
+ l — (For—2 AR) (Ko+D) +R (E +G) D | aÃ 
-+ |— Fo -2 AR) E+G)-R (xo+D) D as ooo hj os Tie (t—t,) | + 
2 32-+(B+)?] ! 
Mm ; 1 
ci Fo — o D)— = ——— — o|— 
+ liam (Fo/—2 AR) (Xo 4 D) R (E+G)D | CT + | — (Fu —2AR) (E+G) + 
po=3(t—t) Cá 
euttmo] B>— me a do sen [(B )(t—to)] + 
B — q 
es ainisiincs Fo, — AR o D e R E CRE ADA a o rea 
+ |[—5 Gu 2AR) Go + D) rola a * | 7 fu-2AR (E+G) — 
= Ra ES ti 
pe BD 4d | e | cos L(B D(t— to)] + 


ae E — 2 + 5?) 1 — 4)Buw 
E) — 2 AR) (D — (Fu— 2 AR)D(E+G) ——————— 
em L + x) Cor aRe Ta UM | E B! + ar 
4 gi 20 (Bº — q?) 
R D)D--——————— — R(EL+G)D--—-—-—--———— — 
+ R (xo + D) (BP — mt)? + 4 3º Bi + R (E+G) (B!— o)? + 4 3 Bº + 


E) [FuA-RD'—-RAS (t—to) + Fo (Fu —2AR)Dcos [2B (t—-t)]—1] + 


To (— FuA-RD'+ RA?) sen [2B (t—to)] + [RG 


-R [(xo + D)'— (E+ G)* | 


—23(t— to) 
— |, € sen (2 (t—to)) + 
4 mo” 


+ [-R60+D) (E R|Go+DE-E+ Ga! 
é] 


* =23)(t=t, mm ER TE me 
: E | A an(t-ml-1] +(E+ t+ 55] e ê toa] 
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(3 em Uy, ser Pi [Gi sen [IB + 0) (t—to)] + Cs cos [(B + 0) (t— ty)] + C: sen [(B—m) (t — to)] 


+ C, cos [(B — nm) (t— 91] + Gs + C (t—to) + C;lcos[2B (t—to)]— 1] + 


25 (t — 


+ Cs sen [2 B (t—to)] + Cs é o “ emulis (eai é 


— 25(t—to)  —23(t-ts) 
+ Co | e “ cos [29 (t—tj—a|+ Cu | e dói cf (37) 


Por substituição dos valores de (36) na expressão (34) verifica-se que 


C = O (38) 
CG + CG=—C (39) 
o que era de esperar, pois equivale a ser o sistema estável quando t—to » co (38) e a ser 
U-—U, = 0 quando t—to==0 (39), Virá finalmente 
u-U,—e t-t)[-(CG+C) sen [(B+0) (t-t)] + Cr cos [(B +) (t—t)] + 
+ Ca sen [(B— mw) (t—to)] + Ci cos [(B— o) (t— to) ] | + C; + 
+ €; [cos [2B(t-to)])— 1]+ Cs sen [2B (t—to)] + 


— 23 (t—t,) — 2 BE 


+ Ge sen (2» (t—to)] + Cwe a [cos[20(t—-t)]-1]) + 


+ Ci [PY eme: Ti) (40) 


Poderemos analisar melhor a expressão (40) determinando valores por saltos Àt que sejam 


a e forçada ça. - Portanto 


fo 


múltiplos simultaneamente dos períodos de oscilação própria ! 


2% 2% 
St=n— =»; — 
B f) 
2 € 2 = 
= (1; — 15) = (n,+n,) = 
( 1 2 pad | 2 Bto 
2 7 
=2n—— = 2h; sa 
2B 200) 


em que ni e ns são inteiros, e portanto ny — ns, ni + n2, 2n1 e 2n2 também o são (*). 


Fazendo então 
7 f == t, = nÃt 
vira 
— Sen. At ZE n.At 


U—Uo = ê CC de Gy le -T Co + Cal] (ui) 


t— t=n8t 


B ' é dá a 
(*) Se E for irracional, é evidente que não haverá dois inteiros que satisfaçam a relação, mas podemos supor 


a sua existência com uma aproximação tão grande quanto quisermos. 
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Em regime permanente (n +09) 
U-—U = €C;— Co— Cu = àU (42) 


AU será positivo, nulo ou negativo consoante a energia potencial inicial é maior, igual ou 
menor que a energia potencial em regime permanente. Se o sistema parte do repouso U, = 0 
e SU >o. Se for perturbado de início por acções do tipo Dirac ou «Step» U, 0 e AUTO. 


A expressão (4) poderá ser parcialmente demonstrada se tomarmos como x a velocidade do sistema 
em regime permanente (28): 


t 
vu, = | [Fo sen [B(t—t)] — Rx] x dt = 
ts 


(er JR (B'— 0,2)? 


m 


“2[(B-o)prasprl 


(Bº — 1,2) r , 
— Bsen [2B (t—t,)] + a, LeosiZB(t-t,)] — 1) | (43) 


e nos mesmos intervalos considerados 
AU = 0 


Note-se que só se demonstra com as expressões (41), (42) e (43) qual o balanço energético 
nos momentos correspondentes a t—t, == n. At sem atender às flutuações que sabemos existir 
nos momentos intermédios. 

A simples observação da expressão (40) garante estabilidade. Com efeito basta verificar que 
na expressão (40) [U — U,] é limitado superiormente. 


1.5 CONCLUSÃO 


Como se verifica, o problema da vibração de sistemas lineares com um grau de liberdade foi 
resolvido por integração directa da equação diferencial. Não houve a preocupação de definir carac- 
terísticas e estudar atributos que, alongando demasiadamente este trabalho, em nada o enriquece- 
riam, uma vez que são correntemente apresentados nas obras da especialidade. 

Deixamos para um trabalho posterior a aplicação do método agora proposto, a sistemas com 
mais de um grau de liberdade. 
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RESUMÉ 


On presente les soluctions usuelles des systêmes vibratoires avec un degré de liberté et 
de caractéristiques lineaires, en trouvant immediatement avec toute generalité des resultats, 
pour intégration par aproximations successives de l'équation ditferentielle correspondante. 
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Estabilidade e repartição dos escoamentos em superfícies 
de aquecimento com circulação forçada 


O — Introdução 


O conhecimento da estabilidade do escoa- 
mento numa superfície de aquecimento é pro- 
blema importante na sua realização industrial e 
assume especial relevo em caldeiras e reactores 
nucleares especialmente nas superfícies evapo- 
ratórias. 

As superfícies são usualmente realizadas por 
circuitos tubulares em paralelo, e o problema 
tem sido extensamente tratado sobretudo quanto 
à estabilidade do escoamento, usando modelos 
físicos que permitem um estudo analítico simples. 
Este é essencialmente de valor metodológico, 
porque o cálculo usual é feito numericamente, e 
por troços. Usando-se frequentemente o com- 
putador electrónico o tratamento usual do pro- 
blema, não só não considera as potencialidades 
deste novo processo de cálculo, como descura 
em geral o problema da repartição dos escoa- 
mentos, entre os vários circuitos individuais 
limitando-se a determinar a estabilidade dessa 
repartição. 

No presente trabalho procura-se considerar 
ambos os aspectos, no âmbito duma teoria não 
vinculada a um modelo físico particular, o 
que permitirá um cálculo prático tanto melhor 
quanto mais próximo da realidade for o modelo 
escolhido, e independentemente de o mesmo se 
traduzir ou não por expressão analítica simples. 

Para tal introduzem-se coeficientes, calculáveis 
por métodos numéricos, que caracterizam o efeito 
da variação dos factores importantes, 

Em trabalho posterior será exemplificada a 
sua aplicação prática, e dar-se-á conta dos mo- 
delos usualmente propostos o que permitirá 
testar a melhoria que se apresenta. 


1 — Proposição do método 


O escoamento em pressão num circuito onde 
se verifica o aquecimento dum fluido que pode 
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sofrer mudança de fase pode em geral ser ex- 
presso sob forma diferencial em função das 
variáveis definidoras do seu conteúdo energético 
específico e das coordenadas espaciais. A intro- 
dução das condições iniciais e aos limites per- 
mite formalmente a descrição da evolução do 
sistema. 

No presente caso interessam-nos apenas os 
regimes estacionários. De igual modo é suficiente 
uma descrição segundo uma linha média do 
escoamento. Introduzindo valores médios no 
espaço (e no tempo), admitiremos que a inte- 
gração das equações conduzirá ao conhecimento 
da perda de carga AP entre as extremidades do 
tubo, em função das restantes variáveis defini- 
doras do estado do sistema, de entre as quais se 
individualiza o caudal mássico G. 

Admitimos pois que se obtém: 


AP=7(6,x1...Xn) 
ou 
f (AP, G, %t 0» yo) = 0 


que define univocamente áP, 

Dado o seu carácter físico admite-se ainda 
que 2(G,xt...xn) ==0 e suas derivadas são 
contínuas, salvo num conjunto numerável de 
pontos de descontinuidade removível, de forma 
que a função pode ser desenvolvida em série 
de Taylor. 

Além disso AP existe, tal como todas as va- 
riáveis no domínio real, e verifica-se que: 


AP(G=0)=0; GD>0o 4P>Do 

1.1 — A função f(AP,G,x1...xn)=0 repre- 
senta uma superfície num espaço afim a n + 2 
dimensões e cada ponto desta superfície define 
completamente o escoamento num circuito. O 
seu conhecimento permitirá resolver todos os 
problemas relativos à associação de vários circui - 
tos, conhecidos um certo número de variáveis 


que individualmente os caracterizam. 
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Designando por ajx os valores da variável xk 
impostos ao circuito j este passa a ser definido 
por 

& (AP, Gy xt, 


reduzindo-se deste modo a descrição a um es- 
paço n42—k dimensões. 

A cada parametrização corresponde uma re- 
dução do espaço de descrição que se reduz a 


[é 


zero quando o escoamento é completamente 


definido. 


. Xjsn-k) == 0 


1.2 — Estabilidade de escoamento num circuito indi- 
vidual 


Sejam xi ==aj...Xn==an Os parâmetros que 
definem o circuito e P*==P* (G) a caracteris- 
tica da bomba que o alimenta. Seja P* a pres- 
são à saída e AP*=P*- Pr=YW(G). 

Os possíveis pontos de equilíbrio do sistema 
são dados por: 


AP = y (G) (1) 
F(G, AP,a...a)=0 (2) 
AP = A Px (3) 


Por hipótese (2) define univocamente AP, 
pelo que eliminando-o entre 1, 2, 3, resulta 


P(G, ar...an) = 4(G) 


que é uma equação em G que terá em geral 
mais de uma raiz. 

Como as variáveis só se definem no domínio 
reale G> O só têm realidade as soluções reais 
e positivas. Estas fornecem as condições neces- 
sárias de equilíbrio do sistema. 

Sejam Gi ... Gr essas raízes. 

Para que nesse ponto o equilíbrio exista e seja 
estável é necessário que para uma perturba- 
ção 3G seja 


(4) AG [9 (G + 4G) — W(G + 4G)]>0 


o que equivale a exprimir que no ponto consi- 
derado a variação de energia mecânica fornecida 
ao fluido pela bomba é menor do que a variação 
da energia dissipada pelo fluido. 

A condição (4) exprime-se fácilmente desen- 
volvendo 7 e y em série de Taylor: 


(5) AG[9(G + 4G) — (Gr +AG)] = 


J ê 2, 
as = ra Rr e my AG' +... 
dG dG |Icgr 211 06º oG* Gr 
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e conclui-se que para que o ponto de equilíbrio 
exista e seja estável é suficiente que o primeiro 
termo da série não nulo seja de ordem ímpar e 
positivo. 

Práticamente, exige se nas aplicações que a 
estabilidade seja determinada pelo primeiro 
termo. 


1.2.1 — A condição de estabilidade determina 
o número de pontos de equilíbrio estável que 
poderão ser em número n>1.Sen>1 eo 
sistema estiver fora dum ponto de equilíbrio o 
algoritmo anterior não permite saber qual o 
ponto estável para que o regime tende. Esta 
indeterminação não é geralmente aceitável em 
condições de funcionamento real de uma su- 
perfície de aquecimento. 

Existindo soluções do sistema de equações, 
a condição suficiente para que n=1 é 
(6) go 

9G JG 
em todo o domínio. 

Praticamente exige-se que (6) se verifique em 
todo o domínio de funcionamento de uma super- 
fície de aquecimento. Subentende-se por isso que 
quando se afirma genêricamente a estabilidade 
se realiza esta condição. 

A relação (6), põe em evidência a característica 
da bomba. Em particular interessa notar que 


se é finito o escoamento é sempre estável 


se -— == — co, sempre instável se Ed = + 00, 
dG dG 

A existência dos dois casos simultâneamente 

não tem realidade prática, porque a caracteris- 

tica duma bomba real não tem pontos de des- 

continuidade infinita. 

O caso mais favorável é praticamente reali- 
zado por uma bomba puramente volumétrica, 
que é aliás a solução frequente quando uma 
bomba alimenta um só tubo. 


É frequente na literatura considerar apenas o 
O 
caso de ser a = 0. Nestas condições, para 


que (6) se verifique e o escoamento seja sempre 
estável é apenas necessário que AP seja mono- 


tônicamente crescente com G. O caso em que 
fosse num ponto =: O é um simples caso 


particular de (5). 


1.2.2 — Na análise anterior definiram-se as 
condições para que um ponto 4Ps, Go, at... an, 
solução do sistema de equações 1, 2, 3, fosse de 
equilíbrio estável. Tem porém interesse conside- 
rar a influência no deslocamento desse ponto da 
variação de um dos parâmetros. Este aspecto 
pode sempre considerar-se repetindo a análise 
anterior para o novo valor do parâmetro. Tal 
envolve porém, praticamente, bastante labor, 
dado que as soluções são geralmente obtidas 
por análise numérica. O problema resolve-se 
mais directamente, especialmente quando as va- 
riações do parâmetro são pequenas, usando a 
teoria das funções implícitas na região de esta- 
bilidade que contém o ponto de partida consi- 
derando o parâmetro como uma variável. 

Suponhamos pois que o parâmetro que pode 
variar é ar e façamos ar =x, sendo (1Po, Go, 
at... Xo--- an) o ponto de equilíbrio inicial. 

O sistema de equações 1, 2, 3, escreve-se, ex- 
plicitando x 

AP==7 (G,x,at...ar—s,ar+1, an) 
AP* = (G) 
AP = AP? 


ou eliminando áP* 

(7) AP — 4 (G)=H (AP, G)=0 

(8) AP —- 26, )=H(Gx)=0 

e o sistema define implicitamente G=G (x) e 
AP = AP (x) pois por hipótese (1) 


d (Fi, fi) 
o — — =D fo 
9) o (G AP) iá 
No novo ponto de equilíbrio será 
dáP 1 /d'AP 
10) AP=AP, —— | Axt.—[ —— jAxte.... 
5 i +( dx my a dx? Ja + 
? 9 
(11) G= Go (SE) a E = Já 
dx /xo 211 dx? /xo 


e como supomos Àx pequeno, podemos despre- 
zar os termos de 2.2 ordem. 


Teremos: 
0 (Ei, ti d» dy 
dAP d (CG, d 
ud SEE (6, x). o dx OG 
dx D oo dy 
dG dG 


(1) Considera-se tal como no que se segue que o 1.º 
termo não nulo em 5, é o primeiro. Como o ponto é de 
equilíbrio estável será pois D > 0, nesse ponto. Além 
disso, supondo-se que a perturbação não afecta a estabi- 
lidade será D > 0 em todo o domínio G, AP, x,ay... an 
de funcionamento, 


4 


OF, fi) sd 

dG dd (x, AP) = dx 
(13) dx o D E dy dy 
dG OG 


e o efeito duma perturbação será tanto menor 
quanto maior 


du 
p= dy 
dG dG 
e quanto menor 
dy dz 
= e ———— — 
9G dx 


que interessa exprimir adimensionalmente em 
ordem a definir um critério comparativo. 
Deste modo definimos: 


fo suit dy ) É LA = S coeficiente de estabi- 
dG Ne SP lidade do escoamento 


, 


a ae = Se — coeficiente de estabilidade 
oG AP do circuito 
so kd = Sp — coeficiente de instabilidade 
AP dG da bomba 
R Mu Sr coeficiente de perturbação 
AP dx 


relativa ao parâmetro r 


Para que um escoamento seja estável é con- 
dição suficiente que S>0 ou Sc >5. 

E o efeito duma perturbação é tanto menor 
quanto maior S e menor Sp e Sr. 


1.2.2.1 — É frequente na literatura [1] caracte- 
rizar um escoamento pelo coeficiente Sc, o que 
corresponde óbviamente a um aspecto restrito. 

Se Sp = 0, tratando-se dum escoamento iso- 
térmico, e para valores muito elevados do número 
de Reynolds é aproximadamente 


AP=k G* 


donde == 2 


E é recomendado em /2] que em qualquer cir- 
cuito aquecido, deve ser S>>1 para que a di- 
ferença máxima na repartição de escoamentos 
havendo circuitos em paralelo ainda seja acei- 
tável. 

Cita-se o valor como regra prática para tubos 
economizadores-vaporizadores de caldeiras, a ti- 
tulo exemplificativo, dado o carácter restrito da 
conclusão e das condições em que é válida como 
se verá quando se tratar de circuitos em paralelo. 


TRONICA 
29 


A conclusão geral é que S >> 0 para que o 
escoamento seja estável, no ponto de funciona- 
mento dado. Se Sp =-0 será em geral S<ã2, 
num circuito aquecido. 


1.3 — Repartição dos escoamentos em circuitos em para- 
lelo quando os colectores não introduzem perda 
de carga 


Seja de novo Fi=-AP*—y (G) e aj; os pará- 
metros que definem os j(j=1...m) circui- 
tos individuais. 

Os pontos que forem de equilíbrio são solu- 
ções de 


(14) E =a4P—-y(G)=0 
m 


(15) Fs=G -—-2G=0 
1 


(16) f =4Pj—o(Gi;,ais) 


| 
fmd 
3 


e tal como anteriormente só podem ser pontos 
de equilíbrio de raizes reais e positivas, que por 
hipótese supomos existirem. A estas podem cor- 
responder pontos de equilíbrio estável ou ins- 
tável, 

Suponhamos que para cada um dos pontos 
(G, Gi, àP) solução de 15 e 16 com ÁPi=4AP é 


9 (Fa, fr... fm) 


am à (O, Sica Ga) 


então o sistema define G = 4 (AP) que é a 
equação dum circuito único equivalente aos cir- 
cuitos em paralelo. 

Calculando o jacobiano (18) atendendo a 14, 
15, 16 e 17 verifica-se que 


7 (Fs r fi; ... Em). mea 
J (G, Gi Ea Gm) 


d Dj 


d G; 


= 3 


para que (18) se verifique é necessário que 


) 
E om 
O G; < 


para todos os j (j="1...m) 


O Tr 


Mas se para alguns dos circuitos é 


<o0 


(r <m) tal significa que nos circuitos indi- 
viduais r, o ponto é de equilíbrio instável, donde 
o sê-lo para o circuito global em paralelo. 
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Para que o escoamento seja estável é pois 
necessário que 


(19) 
9G; 
o que reduz o problema da estabilidade no para- 
lelo ao dos circuitos individuais alimentados por 
uma bomba de características 4 ==const. Deste 
facto o interesse de conhecer Sjc, para cada 
tubo. 
Verificando-se (18) é 


A 
(20) gg nes! 
dG q À 
j 09 
Gi 


e verificando-se (19) o escoamento global é está- 
vel, e o coeficiente de estabilidade do paralelo, 
é o do circuito global equivalente: 


TR 
y 0d Dj dG AP 
j dG 


1.3.1 — Verificando-se (18) e (19), existe um 
circuito global equivalente que é estável. Neste 
caso, o sistema de equações 14, 15, 16, 17 define 
a repartição dos escoamentos no ponto de equi- 
líbrio. 

Interessa considerar de novo o efeito duma 
perturbação num parâmetro dum dos circuitos 
individuais naquela repartição, além do seu efeito 
no escoamento global mantendo-se ainda estável 
o escoamento (1). O problema trata-se como em 
1.2, mas repete-se por forma a incluir aquele 
como um caso particular das expressões obtidas. 

Seja o sistema de equações, cujas soluções con- 
têm todos os pontos de equilíbrio, que se supõe 
existirem e serem estáveis: 


(21) B=AP—y4(G)=0 
m 

(22) E=G —26G=0 
l 


(23) fi = AP — q; (G;)) = 0 (is lc 
(24) f=AP-»[(G,x)=0 
onde se designa por x o parâmetro a;;, do cir- 
cuito |, que sofreu uma variação. 
Õ ponto (A Po, Go, Giro ee. Gno, Xo) é por 
hipótese de equilíbrio. 
d (F . 
Como (Fu, Fa, fr... Em) D>0 
o (AP, G, Gi... Gn) 
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O reactor de água fervente 


General Electric 
tornar-se-d o seu 


SR o des sao dedo 


3 BILIÕES KILOWATTS-HORA 


Durante Outubro 1963 cinco centrais com reactores 
de água fervente General Electric, produziram mais 
de 3 biliões de kWh, a maior parte nos últimos 
30 meses. Em 1968, pelo menos mais 6 reactores G-E 
do mesmo tipo estarão em actividade produzindo mais 
do que 16 biliões de kWh. 

Falando de experiência — experiência operacio- 
nal — mostramos as vantagens da simplicidade do 
desenho do reactor G-E de água fervente, com bai- 
xos custos operacionais, fornecimento e manutenção 
fáceis. 


mi investimento 


O reactor simples, de ciclo directo precisa de menos 
material e menos peças auxiliares. Não há permu- 
tadores de calor que possam complicar a manutenção 
do sistema. O vapor é fornecido directamente para 
a turbina General Electric a 1000 psi e 545º F. 

O edifício do reactor é acessível mesmo durante 
o seu funcionamento. Para simplificar a manutenção 
do reactor e o seu abastecimento com combustível, 
a G-E usa um sistema controlado mecânicamente, 
que torna mais eficiente a utilização de combustível. 
Segundo a opinião dos operadores da Central 
Nuclear de Dresden (220,000 kWe) esta central é, 
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dido es 
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] o w a 
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E des 


até ao presente, a melhor unidade de produção de 
vapor. 

Quando a turbina precisa de ser inspeccionada, a 
experiência mostra que não há qualquer dificuldade. 
Logo que a turbina arrefeceu suficientemente para se 
poder trabalhar nela, pode-se fazê-lo sem qualquer 
receio pois a sua rádio-actividade é nula. 

Os reactores de água fervente General Electric são 
presentemente os mais económicos produtores de 
energia de origem nuclear. A General Electric está 
oferecendo reactores de água fervente com qualquer 
capacidade que se pretenda, a preços fixos, e com 
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garantia de fabrico. A General Electric gostaria de 
ter uma oportunidade de discutir este assunto com 
os interessados. Contacte a General Electric Portu- 
guesa, Rua do Norte, 5, Lisboa-l, ou escreva para: 
International General Electric Company, Dept. 30-32 TC, 
159 Madison Avenve, New York, NY 10016, U.S.A, 
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por hipótese, o sistema (21), (22), (23) e (24) 
define implicitamente AP, G, G;, em função de 
(x) na vizinhança do ponto. 

Quando xo passa a xo + Ax o novo ponto 
de equilíbrio será (1) 


Au 
E 
"q 


AP= ap +(SO) Axo 


dx 

G = Go +() dutos 
dx 

G; = Go + (4) xt. 
dx 


em que de novo se consideram apenas os pri- 
meiros termos por se supor Àx pequeno. 


Será d (Fu, Fa, fi... Em) 
daP  d(x,G, Gr... Gn) 
dx D 
d (Ei, Es, fy ... Em) 
dG | Od(AP,x,Gi...Gm) 
dx D 
o (Fr, Fs, £4 ... fu) 
dG —  o(AP,G,Gi... x... Gm) 
dx D 


Atendendo a 21, 22, 23 e 24, vem: 


es o (E, Es, fi ... Em) Sae 
d (AP, G,G4 ... Gm) 


m | dp m 17 
1 dG 1 07] 
dG; 


e será: 


(66 G Go) 
— dn dy n dj 
dx OG 1 9G; 


d = m 
Oto tm) Um dA 4 mf 
à (AP,x, Gi... Gm) dx À OG; 


(j=1..m m£)) 


e resulta: 
dg 
dáP Ox 1 
de dg 'fãa ção 
9G; ! d9j NE dy 


1 
dj dy 
dG dx dG 
dx 031 moto EK 
oG; : “09 E dy 
dG; dG 
d 04 
dGk = dx 
dx dp Ok 
0G: dGk 
x — E 4a | m =! 
sã 1 
10%; OU 
E dE. 
O 
du | DE 
dG dx E 41 
dx dy | dg, (: 1 ag) 
dG: | 0G; 1 0% dy 
E oG; dG 


1.3.2 — As expressões obtidas generalizam as 
expressões para um só circuito e permitem defi- 
nir, analogamente, mas com toda a generalidade 
os coeficientes de estabilidade e de perturbação. 

Neste caso, define-se mais o coeficiente de 
perturbação no escoamento do circuito j por 
perturbação do parâmetro r no circuito |, como 


o q! 
5 x X* | (11) 
é dq : d Dj G 
G: 06; 


analogamente se definem os restantes coeficientes. 

Facilmente se verifica a influência da caracte- 
rística da bomba, responsável pela transmissão 
duma perturbação num circuito aos restantes, 
donde o aspecto restrito das conclusões que 
reduzem a análise ao de um circuito isolado. 

Como o demonstra a análise anterior, só no 
caso duma característica Y (G) = const o efeito 
duma perturbação se reduz ao circuito individual 
onde se verifica. 

A análise clássica junta a esta simplificação a 
de que os circuitos são praticamente iguais, e 
desconhece o efeito de perturbações que existem 
sempre num escoamento real. Destas assinála-se 
como particularmente importante o que se refere 
ao aquecimento dos circuitos individuais. 


TRONIGA 
31 


2 — Conclusão 


A utilização de expressões analíticas no cál- 
culo da estabilidade e da repartição dos escoa- 
mentos conduziu clâssicamente à adopção de 
modelos simplificados, dos quais citaremos o 
mais usado, em próximo artigo. 

Acentua-se porém, que a análise dum caso real 
não implica para a sua solução o manejo de 
expressões analíticas complexas, como se poderia 
induzir duma teoria que pressupõe a existência 
das várias funções utilizadas. 

O cálculo prático é feito numericamente cal- 
culando AP, em função de Gi, e fazendo o 
cálculo por troços. Deste modo e apenas na vizi- 
nhança do ponto que corresponde ao regime de 
interesse se determina o circuito equivalente, 
No processo usual, determina-se neste ponto Sc, 
e procura-se que seja superior a um valor mi- 
nimo. Se o não for procura-se a sua correcção, 
de modo que será exemplificado posteriormente. 

O cálculo assim feito, que se baseia em mode- 
los simples, torna extremamente complexa e tra- 
balhosa a exploração duma variação nos pará- 
metros, porque na forma usual de proceder 
equivale a uma nova determinação do ponto de 
equilíbrio, ou seja a uma repetição sucessiva do 
mesmo problema. 

O método proposto permite resolver com faci- 
lidade o problema das pequenas perturbações em 
torno de qualquer ponto de equilíbrio, que é 


muitas vezes o caso importante, e o formulário 
aplica-se a toda a região onde existe equilíbrio, 
sem que as suas conclusões estejam vinculadas 
a um modelo físico particular. Para tal, como é 
apenas necessário o cálculo de primeiras deriva- 
das, o mesmo efectua-se numeéricamente, e o es- 
tudo faz-se com toda a generalidade a partir 
dum número restrito de parâmetros com as rela- 
ções encontradas. 

Dispondo dum computador, o aspecto fastidioso 
do cálculo é eliminado, e torna-se possível pre- 
ver a influência de cada perturbação, e com a 
maior generalidade. 

O método responde deste modo às possibili- 
dades dos novos meios do cálculo numérico, e 
integra imediatamente qualquer melhoria no 
conhecimento da determinação da perda de carga, 
sem maior complicação que a que resulta dum 
calculo numérico por troços. 
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RESUMO 


Admitindo propriedades muito gerais quanto à dependência da perda de carga 
das variáveis definidoras do escoamento num circuito aquecido, estabelecem-se as 
condições de estabilidade num circuito individual e em circuitos em paralelo. Linea- 
rizando a função num domínio de estabilidade determina-se a influência da perturbação 
num dos parâmetros em torno dum ponto de equilíbrio. Introduzem-se critérios adi- 
mensionais de estabilidade do circuito, da bomba e da perturbação. 


SUMMARY 


With some general assumptions about functional dependence of pressure drop of 
a heated tube on the flow variables the conditions for stability of flow on individual 
and parallel tubes are deduced. Adimensional critéria are introduced for stability and 


perturbation. 


TLONICA 
32 


PoclaimnT.Y.45 


Escavadora-Retro Escavadora-Grua 


e POTÊNCIA AUMENTADA 

e ROTAÇÃO TOTAL E CONTÍNUA 
e ESTABILIDADE 100º 

e 20 BALDES DIFERENTES 


REPRESENTANTES EXCLUSIVOS 

Sociedade de Mecanização Industrial e Agrícola, S.A R.L. 
Avenida Padre Manuel da Nóbrega, 8-B-1.º — LISBOA 
Telefone: 724053 4/5 


TECNICA XXVII 


impresa Ultramarina de Sondagens E Fundações, |." 


Associada: Empresa de Sondagens e Fundações 
TEIXEIRA DUARTE, LDA. 


SONDAGENS geológicas e geolécnicas — CAPTAÇÕES de águas 

subterrâneas — REBAIXAMENTOS do nivel freático — CONSOLIDA- 

ÇÕES E ESTABILIZAÇÃO do solo — INJECÇÕES de cimento e 

outros produtos — INFRAESIRUTURAS de barragens, túneis e pontes 
mes FUNDAÇÕES de todos os tipos 


Única Empresa da especialidade com sede no Ultramar Português 


Sede-LUANDA 
Rua Serpa Pinto, 60 LISBOA L. MARQUES =C. P. 932 
C. P. 2178 — Tel. 4226  P.da Figueira, 18, 3.º Esq. BEIRA — C. P. 1002 


Companhia Portuguesa de Trata || --;MPERMEÁBILIS” 


S.A. R.L. 
ARMADURAS PRE-FABRICADAS MALHASOL O IMPERMEABILIZANTE OFICIALMENTE 
ACONSELHADO EM TODAS AS ARGA- 


| MASSAS DE CIMENTO 
Prédios, empenas, caves. — Depósitos para água, vinhos 
e aguardentes. — Facultam-se certificados de ensaio 
dos Laboratório Nacional de Engenharia Civil e da 

Junta Nacionai do Vinho 


T a - e s , z 
Fensão de segurança 3.000 Kg em Distribuldor geral para todo o Império Português: 


AÇO Bi ANSELMO DE MATOS 
Av. Almirante Reis, 179, r/c. - Tel. 46439-LISBOA - 1 


Revendedores em Lisbos: 


João Pereira Vareiro & Filhos, Lda. 
Rua Augusto José Vieira, 12 - Tel. 847480 


“Manuel Pereira Matias, Lda. 


Tensão de segurança 4.000 Kg'cm? Rua de Pedrouços, 105-A - Tel. 6r 1153 
AÇOS DE ALTA RESISTÊNCIA PARA BETÃO ARMADO Manufacturas de Cimento Fortex, Lda. 
de Toqieseregata tésaioa é económica a que conduzem são Etna der Bonfim Jeje io 
Para qualquer obra em betão armado comsulte: Sequeira & Santos Lda 
, . 


DIAL — GABINETE TÉCNICO 


R. DO INSTITUTO INDUSTRIAL, 18-1,º DT.º —- LISBOA 
TELEFONES 6712245 E 673785 ? 


R. José Joaquim Marques, 113- Tel. 0261 05 - MONTIJO 


TÉCNICA - XXVII 


C. D. U. 539.443: 539.43 


FLEXÃO SIMPLES 


Contribuição para a apresentação do cálculo em relação à rotura (Fase |l1) 


«Ampliemos os nossos conhecimen- 
tos sem a ideia fixa de encontrar 
cálculos tipo receita.» 


Saliger 
1 — INTRODUÇÃO 


1.1 — Entre nós, infelizmente, não está sendo 
utilizado, em relação à Fase III, o processo de 
determinação da resistência ou seja, a verifica- 
ção de secções e cálculo da armadura; tal facto 
afigura-se-nos de lamentar, em face da intensa 
evolução sofrida pelo estudo das peças de betão 
armado, evolução que teórica e experimental- 
mente comprovou a falta de justeza do cálculo 
em Fase II, ou seja, o cálculo clássico, principal- 
mente nos casos em que o esforço é maior e 
portanto quando o «rigor» e economia se torna- 
riam mais necessários. Deve frizar-se que esta 
vantagem do cálculo em relação à rotura sobre 
o cálculo clássico é tanto maior quanto maior 
for o esforço; assim, se nas peças pouco esfor- 
çadas essa vantagem mal se notaria, nas muito 
esforçadas, seria surpreendente a economia veri- 
ficada (vejam-se exemplos no final deste traba- 
lho). Se por outro lado quisermos atender, de 
modo objectivo e particular, à qualidade dos 
materiais a empregar, o «cálculo à rotura», tor- 
na-se muito mais realístico e seguro que o clás- 
sico. Os métodos clássicos, assentando em bases 
elásticas e portanto em diagramas tensões-defor- 
mações lineares, «ignoram» as possibilidades 
plásticas dos materiais, não permitindo assim 
uma exploração completa e perfeitamente viável 
das suas qualidades. 

Como justificação do esquecimento a que 
parece ter sido votado o cálculo à rotura, é cos- 
tume responsabilizar-se o Regulamento do Betão 
Armado, dizendo-se que ele não o permite, o 
que aliás estaria de certo modo em contradição 
com o Regulamento de Segurança das Constru- 


por G. FARIA DE VASCONCELLOS 


(ENG. CIVIL) 


«O objectivo do cálculo estático é mais determinar ou prever 
a segurança duma estrutura relativamente à rotura, de que 
determinar as tensões nela causadas pela solicitação prevista,» 


E. Mórsch (1909) 


ções Contra os Sismos visto nele ele estar pre- 
visto. Sejamos justos quanto ao R.B.A., pois que 
ele simplesmente omite, como aliás seria de pre- 
ver, se atendermos à data da sua promulgação e 
omitir, neste caso, não é proibir. O nosso R.B.A. 
diz textualmente : «os cálculos de resistência não 
podem ser feitos por métodos empíricos mas 
somente por métodos científicos baseados em 
dados experimentais»; ora o cálculo em relação 
à rotura tem precisamente uma base científica e 
experimental mais desenvolvida e actualizada 
que o cálculo clássico, conduzindo a resultados 
praticamente coincidentes com os obtidos expe- 
rimentalmente o que, nos casos mais delicados, 
não se verificava com o cálculo em Fase II. O 
espírito dum regulamento de carácter técnico, 
com certeza não é, nem pode ser, O contrariar ou 
entravar o desenvolvimento e aperfeiçoamento 
dos processos de cálculo (que é o caso) ou de 
construção, mas sim, «ajudar» os técnicos e evitar 
que se ultrapassem ou atinjam valores tidos como 
arriscados ou perigosos, fornecendo implicita- 
mente valores e conclusões que só a investiga- 
ção laboratorial o poderia fazer. Na apreciação 
comum chega a parecer que entra em jogo um 
efeito psicológico: a palavra rotura, que ame- 
dronta e que possivelmente exige assim outros 
betões e outros cuidados, quando afinal o nome 
«cálculo à rotura» é abreviação deste outro 
«cálculo em relação à rotura», o que, por si só, 
não implica o emprego de betões melhores nem pio- 
res dos que usualmente se empregam. 
Entretanto será conveniente esclarecer desde 
já, que, seja qual for o «cálculo em relação à 
rotura», em betão armado, é factor obrigatório 
não tolerar qualquer espécie de plastificação 
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